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Abstract 

The ground state structures of small sili- 
con clusters are determined through exhaus- 
tive tight-binding molecular dynamics sim- 
ulation studies. These simulations revealed 
that Sin is an icosahedron with one missing 
cap, Sii 2 is a complete icosahedron, Si 13 is 
a surface capped icosahedron, Siu is a 4-4-4 
layer structure with two caps, Siis is a 1-5- 
3-5-1 layer structure, and Siig is a partially 
closed cage consisting of five-membered rings. 
The characteristic feature of these clusters is 
that they are all surface. 

Smalley and co-workers discovered that 
chemisorption reactivities of silicon clusters 
vary over three orders of magnitude as a 
function of cluster size. In particular, they 
found that Si 33 , Si 39 , and Si 45 clusters are 
least reactive towards various reagents com- 
pared to their immediate neighbors in size. 
We provide insights into this observed reac- 
tivity pattern through our stuffed fullerene 
model. This structural model consists of 
bulk-like core of five atoms surrounded by 
fullerene-like surface. Reconstruction of the 
ideal fullerene geometry gives rise to four- 
fold coordinated crown atoms and 7r-bonded 
dimer pairs. This model yields unique struc- 
tures for Si33, Sisg, and Si45 clusters with- 
out any dangling bonds and thus explains 
their lowest reactivity towards chemisorption 
of closed shell reagents. We also explain why 



a) these clusters are substantially unreactive 
compared to bulk surfaces and b) dissocia- 
tive chemisorption occurs on bulk surfaces 
while molecular chemisorption occurs on clus- 
ter surfaces. Finally, experiments on Si^X^ 
(X = B, Al, Ga, P, As, A1P, GaAs) are sug- 
gested as a means of verifying the proposed 
model. 



1 Introduction 

Nanometer size particles are the embryonic 
forms of matter whose microscopic study pro- 
vides insights into the evolution of mate- 
rial properties from molecules and surfaces 
to solids |lj, 0, |J. Furthermore, nanoscale 
particles have been shown to exhibit exotic 
optical properties and reactivities quite dif- 
ferent from those in molecules and solids 
0, 0, |, § §, §. For these reasons, theo- 
retical studies on clusters are critical to the 
design and synthesis of advanced materials 
with desired optical, electronic, and chemi- 
cal properties. Such studies are at the in- 
terface of the traditional fields of quantum 
chemistry, solid state chemistry, and statis- 
tical mechanics. Hence, physicists, chemists, 
and material scientists are working individ- 
ually and in teams to unearth the funda- 
mental principles underlying the structure, 
dynamics, and reactivities of these clusters 
|, |, |, |, |, |, 0, 



ljpira. Indeed 
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the diverse group of scientists assembled at 
this Nanomeeting-95 is testimonial to the 



1 



breadth and depth of scientific inquiry into 
this novel state of matter. 

The nanoscale semiconductor clusters have 
potential applications in non-linear optical 
devices, photovoltaic devices, and as photo- 
chemical catalysts. For this reason, intense 
experimental activity exists on the synthe- 
sis, structural characterization, spectroscopy, 
and reactivities of these clusters. On the 
other hand, the theoretical and computa- 
tional efforts on these systems are at a prim- 
itive stage. In particular, due to both exper- 
imental and theoretical limitations, the de- 
tailed knowledge of the role of surface struc- 
ture on the reactivities of these nanoparti- 
cles is lacking. To correct this imbalance, 
we started a systematic theoretical investi- 
gation of the structure-reactivity relationship 
in semiconductor clusters. The results from 
these investigations is the subject of this Ar- 
ticle. 



2 Structures of small sili- 
con clusters 

Since naked silicon clusters are highly reac- 
tive, they are mostly synthesized in a molecu- 
lar beam under high vacuum conditions || § . 
The number density of available clusters is 
so low that diffraction based structural in- 
vestigations are not feasible under these ex- 
perimental conditions. Theoretical calcula- 
tions have now established the structures 
of Si 7v clusters in the N = 2-10 atom size 
These calculations have em- 



range |14, |15 



ployed quantum chemistry molecular orbital 
techniques. While these methods are accu- 
rate for structural determination, computa- 
tionally they are highly cpu intensive and 
hence it is difficult to search all possible clus- 
ter geometries. Consequently, we employed 
the tight-binding molecular dynamics (TB- 
MD) simulations for the determination of the 



ground state structures of Si^ clusters in the 
iV = 11-16 atom size range. The results of 
these simulations are presented below. 

The choice of the tight-binding method for 
the study of the cluster structures is moti- 
vated by its accuracy and computational effi- 
ciency jn| [17|, |18|, [19]]. This all valence elec- 
tron method is equivalent to the extended 
Hiickel method well known in theoretical 
chemistry [pO |. This method has yielded 
structures in excellent agreement with ab ini- 
tio electronic structure calculations for both 



carbon and silicon clusters [16], [L7|, [18], [19], 
Recently, Menon and Subbaswamy have con- 
structed an accurate tight-binding Hamilto- 
nian for silicon clusters |I8|, Il9||. This Hamil- 



tonian includes Harrison's universal param- 
eters appropriate for the description of bulk 
Si 
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supplemented with two to four addi- 
tional parameters for the description of the 
silicon clusters. These additional parameters 
are derived by fitting to the Si2 bond length 
and vibrational frequency and to the overall 
size-dependent cohesive energy curve of clus- 
ters [P| [Tj|. Most importantly, none of the 
parameters of this Hamiltonian are fit to any 
of the cluster structures. Full details of the 
Hamiltonian and computational methods are 
described elsewhere fjl| [19| . 

As the size of the cluster grows, the num- 
ber of structural isomers increases exponen- 
tially, with the result that searching the com- 
plete configuration space for the global poten- 
tial energy minimum becomes a formidable 
task. However, by combining the tight- 
binding method with the molecular dynamics 
|23[| simulated annealing technique we can ef- 
ficiently search the cluster configuraton space 
and determine the ground state geometry 
2|, We used this TB-MD method in all 



the calculations reported here. 

Figure 1 displays the lowest energy struc- 
tures obtained using this TB-MD technique 
for clusters in the N = 11-16 atom size range. 
The Sin structure consists of two tetragons in 
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the anti-prism geometry and three caps. Two 
of these caps are attached to the opposite 
faces of the two tetragons, while the third one 
is attached to the edge of the top tetragon. 
This structure may also be described as an 
icosahedron with one missing atom. Replac- 
ing this missing atom would give rise to the 
icosahedral cage structure for Sii 2 - Such a 
spherical cage structure has not been pre- 
dicted or observed for any of the 12-atom el- 
emental clusters. Adding a face cap to Sii2 
gives the lowest energy Sii 3 structure. An al- 
ternative structure, derived by placing a Si 
atom inside the cage of Sii 2 , is a high energy 
local minimum. 

Adding a face cap to the Sii3 structure does 
not yield the lowest energy Sii 4 structure. In- 
stead, Sii4 assumes a layer structure consist- 
ing of three planes of four atoms each and 
two adjacent face caps. By suitably rotat- 
ing this structure, we may also describe it 
as a pentagon sandwich (or prism) with two 
caps each at the top and the bottom. The 
pentagon prism is somewhat distorted and 
displaced. We also considered a bi-capped 
hexagonal anti-prism as a candidate for the 
ground state of Sii 4 . However, this struc- 
ture proved to be unstable, indicating that 
six-atom ring structures are still not favored 
in these small clusters. Like Siu cluster, Siis 
assumes a layer structure consisting of 1-5-3- 
5-1 layers. However, Siie is quite unlike any 
of the previous clusters. It is an open cage 
consisting of fused pentagons, reminescent of 
the small fullerenes. 



3 Stuffed Fullerene Model 

Searching the complete configuration space 
for the ground state structures is very dif- 
ficult for large clusters, even with the TB- 
MD method. Consequently, we propose a 
simple model that explains the trends in the 
reactivities of silicon clusters in the 30-50 



atom size regime. We call this the stuffed 



fullerene model [26|. This model consists of 
1) a central atom A, 2) four atoms B sur- 
rounding A in tetrahedral geometry, and 3) 
fullerene surface. The B atoms bind to twelve 
surface atoms, thus rendering the B atoms 
also bulk-like with four-fold coordination and 
tetrahedral geometry. The surface then re- 
laxes from its ideal fullerene geometry, the 
same way the (lxl) bulk surfaces relax. 
This relaxation gives rise to crown atoms and 
dimers (CAD) pattern on the surface. The 
crown atoms are formally three-fold coordi- 
nated and possess one dangling bond each. 
The dimers are also formally three-fold co- 
ordinated, but they eliminate their dangling 
bonds through 7r bonding. The essential fea- 
ture of this construction is that the bulk-like 
core of five atoms (A + B) and fullerene-like 
surface make these structures stable. In prin- 
ciple, this model is applicable to clusters con- 
taining more than twenty atoms. 

Unlike carbon, silicon does not form strong 
delocalized 



7r bonds 
cage structures 



Consequently, 
are 



fullerene 

energetically unfavorable for silicon clusters 
3"H| . Instead, intermediate sized silicon clus- 



ters prefer cr-bonded network structures sim- 
ilar to the diamond structure of bulk sili- 
The fullerene geometry for the sur- 



con 
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face, consisting of interlocking pentagons and 
hexagons, gives special stability to the surface 
atoms p9l. Furthermore, since delocalized n 



bonding is not favorable in silicon, we expect 
the surface atoms to relax from their ideal 
fullerene geometry to allow for dimer forma- 
tion through strong local 7r bonding. Our 
model accounts for all these facts. 

We generate the Si33 structure by stuffing 
the 5-atom core inside the Si28 fullerene cage. 
We orient the 5-atom pyramid in such a way 
that the central atom A, the core atom B, 
and the crown atom C lie on a line. The 
crown atom is at the center of three fused 
pentagons and it is surrounded by three other 
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surface atoms D. The D atoms now relax in- 
wards to form the B-D bond. The same type 
of relaxation motion is necessary to form the 
2x1 reconstruction on the bulk Si(lll) sur- 
face |32|, P3f| . The activation barrier of ~ 0.01 
eV |32], |34| for this relaxation is easily re- 
covered by the formation of the B-D bond, 
whose strength is 2.3 eV/bond |27], |33|. Con- 
sequently, such a relaxation of fullerene sur- 
face is feasible even at 100 K. Finally, the 
remaining surface atoms E readjust to form 
as many dimers as possible. Similar relax- 
ation motion and dimer formation occurs on 
Si(100) surfaces also [^, [J5J. The dimers 
are a-bonded pair of atoms whose dangling 
bonds are saturated through the formation 
of 7r bonds. 

The structures of Si 33 , Si 3g , and Si 45 clus- 
ters thus generated are displayed in Fig. 2. 
These structures reveal that the crown atoms 
are able to form a fourth bond to the core 
atoms B, thus rendering the B atoms for- 
mally five-fold coordinated. The B-C bond 
arises from the back donation of the electrons 
from C to B and it weakens the neighboring 
bonds through electronic repulsion. We do 
not know the strength of this back bond, but 
it is sufficiently strong to eliminate the dan- 
gling bond on the crown atom and make these 
magic number clusters unreactive. 

Chemisorption is sensitive to the local elec- 
tronic structure and chemical bonding on the 
surface. On bulk Si(lll)-(7 x 7) surface, 
chemisorption occurs preferentially on rest 
atom sites and less on crown atom sites 
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Absence of rest atoms in silicon clusters is the 
main reason why these clusters are much less 
reactive than bulk silicon. Furthermore, dis- 
sociative chemisorption of ammonia on bulk 
Si(lll)-(7x7) surface is due to the close prox- 
imity of the rest atom and crown atom sites 
|313fl . Absence of such a configuration on sil- 
icon clusters explains why reagents are un- 
able to dissociate on the cluster surfaces; in- 
stead they chemisorb molecularly. We predict 



that the alloy clusters Si x X y (X = B, Al, Ga, 
P, As, A1P, GaAs) will be highly inert and 
that it may be possible to synthesize these in 
macroscopic quantities. We suggest experi- 
ments on these alloy clusters as a means of 
fully exploring the role of surface structure 
on reactivities. 



4 Summary 

In summary, we determined the structures of 
small silicon clusters through TB-MD simula- 
tions. These simulations revealed that Sin is 
an incomplete icosahedron, Sii2 is a complete 
icosahedron, Sii 3 is a surface capped icosa- 
hedron, Sii4 is a 4-4-4 layer structure, Sii 5 
is a 1-5-3-5-1 layer structure, and Sii 6 is a 
partially closed cage consisting of fused pen- 
tagons, reminescent of small fullerenes. The 
atoms in all these clusters strongly prefer to 
lie on the surface rather than inside. 

We also propose a novel structural model 
to explain the experimental data of Smal- 
ley and co-workers on the chemisorption re- 
activities of silicon clusters towards various 
reagents ||. This model consists of bulk-like 
core of five atoms stuffed inside reconstructed 
fullerene cages. The resulting structures of 
Si 33 , Si 39 , and Si 45 are unique, have maxi- 
mum number of four-fold coordinated atoms, 
minimum number of surface atoms, and zero 
dangling bonds. This model does not yield 
such unique structures for other intermediate 
sized clusters and hence they will have larger 
number of dangling bonds. This explains why 
Si 33 , Si 3 g, and Si45 clusters are least reactive 
towards reagents with closed shell electronic 
structure, such as ammonia, methanol, ethy- 
lene, and water ||. We suggest that exper- 
iments on alloy clusters are needed to fully 
understand the role of surface structure on 
the chemical reactivities of these nanometer 
size particles. 
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Figure 1: The lowest energy structures of Si^v 
(N = 11-16) clusters, determined using the 
tight-binding molecular dynamics simulated 
annealing technique. 



Figure 2: Structures of the least reactive clus- 
ters, Si 33 , Si 39 , and Si 45 , obtained using the 
proposed stuffed fullerene model. These clus- 
ters do not possess any dangling bonds and 
hence are least reactive towards reagents with 
closed shell electronic structure, such as am- 
monia, methanol, ethylene, and water. Rep- 
resentative atoms in different chemical envi- 
ronments are labeled from A to E. 
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